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Abstract 

Motivated by a number of recent experimental studies, we have carried out the microscopic 
calculation of the quasiparticle self-energy and spectral function in a doped graphene when a 
symmetry breaking of the sublattices is occurred. Our systematic study is based on the many- 
body GqW approach that is established on the random phase approximation and on graphene's 
massive Dirac equation continuum model. We report extensive calculations of both the real and 
imaginary parts of the quasiparticle self-energy in the presence of a gap opening. We also present 
results for spectral function, renormalized Fermi velocity and band gap renormalization of massive 
Dirac Fermions over a broad range of electron densities. We further show that the mass generating 
in graphene washes out the plasmaron peak in spectral weight. 
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I. INTRODUCTION 



Graphene is a single atomic layer of crystalline carbon on the honeycomb lattice consists 
of two interpenetrating triangular sub lattices A and B, has opened up a new field for funda- 
mental studies and applications. jl|, 2, [s, Peculiar electronic properties of graphene give 
rise possibility to come over silicon-based electronics limitations. 5|] The single-particle en- 
ergy spectrum in graphene contains two zero-energy at and K~ points of the Brillouin 
zone which are called as valleys or Dirac points. Due to the presence of the two carbon 
atoms per unit cell, the quasiparticle (QP) need to be described by a two component wave 
function. 

The charge carriers in a pristine graphene show linear and isotropic energy dispersion re- 
lation and massless chiral behavior for the energy scales up to 1 eV. Recently, graphene has 
revealed a variety of unusual transport phenomena characteristics of two-dimensional (2D) 
Dirac Fermions such as an anomalous integer quantum Hall effect at room temperature, a 
minimum quantum conductivity, Klein tunneling paradox, weak and anti-localization, an 
absence of Wigner crystallization phase and Shubnikov-deHaas oscillations that exhibit a 



phase shift of vr due to Berry's phase 
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12| One important difference between 



conventional electron gas and Dirac Fermion particle is that the contribution of exchange 
and correlation to the chemical potential is an increasing rather than a decreasing function 
of carrier-density. This property implies that exchange and correlation increase the effective- 
ness of screening, in contrast to the usual case in which exchange and correlation weakens 
screening. This unusual property follows from the difference in sublattice pseudospin chiral- 
ity between the Dirac model's negative energy valence band states and its conduction band 
states. 

The massless Dirac-like carriers in graphene have almost semi-ballistic transport behav- 
ior with small resistance due to the suppression of back-scattering process, and moreover 
graphene is a good thermal conductor. ISj The mobility of carriers in graphene is quite 
high [isl, Q, [17 1 which is much higher than the electron mobility revealed on the semi- 



conductor hetrostructures. 



18 



191 ] On the other hand, by measuring the stiffness of materials 



it is shown that graphene is the strongest material in two-dimension structures. 20|, l21| These 



properties as well as capability to control 



voltage or the chemical doping 22 



23 



24 



of the type and density of charge carriers by gate 



25| make graphene an ideal candidate for superior 
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nano-electronic devices operating at high frequencies. 

Most electronic apphcations are based on the presence of a gap between the valence and 
conduction bands in the conventional semiconductors. The band gap is a measure of the 



threshold voltage and on-off ratio of the field effect transistors (FETs). [26|, |27| Therefore, 
for integrating graphene into semiconductor technology, it is crucial to induce a band gap 
in Dirac points to control the transport of carriers. Consequently, band gap engineering 



in graphene is a hot topic with fundamental and applied significance. [28|] In the literature 
several routes have being proposed and applied to induce and control a gap in graphene. One 
of them is using quantum confined geometries such as quantum dots and nanoribbons. 29|, 
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33| It is shown that the gap values increases by decreasing of nanoribbon width. 



Another alternative way is spin-orbit couplin g wh ose origin is due to both intrinsic spin-orbit 



interactions and the Rashba interaction. 
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371] Another method to generate a gap 



in graphene sheets is an inversion symmetry breakingof the sublattices when the number 



of electrons on A and B atoms are different 
graphene on proper substrates 
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411 or Kekule [42] distortion, e.g. 
48| or adsorb of some molecules 



J44 

52| and CrOg [53| or an alkali-metal sub-monolayer on graphene 



such as water, ammonia 
sheets. 

Recently angle resolved photoemission spectroscopy ( ARPES) experiments on graphene 
epitaxially grown on SiC and ab initio simulations reported a gap opening in the band 
structure of graphene placed on proper substrates, and suggested that interactions between 



the graphene sheet and the substrate leads to symmetry breaking of the A and B sub 



and it consequences to induce a gap in the band structure. Experimenters 24, 



25|, 



43 



attices 



44 



54] 



observed a gap of 260 meV in band structure of the epitaxi al g raphene on SiC substrate 



due to interaction with substrate. In addition, Zhou et al. [2J] found a reversible metal- 
insulator transition and a fine tuning of the carriers from electron to hole by molecular 
doping in gapped graphene. A Density Functional Theory (DFT) calculation confirmed a 
substrate induced symmetry breaking. 55|]. Their results showed a gap in the band spectra 
of graphene about 200 meV which is in agreement with recent experimental observation. 
Their calculation determined that there is a 140 meV on-site energy difference between 
two sublattices. In addition, a band gap is observed in spectra of graphene on Ni(lll) 



substrate 
substrate 



45 



461 ] as well as a gap about of 10 meV in suspended graphene above a graphite 



471 ] due to sublattice symmetry breaking mechanism. Moreover, based on the ah 
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initio calculations, it is suggested that boron nitride substrate induced a gap of 53 meV. [48 1 
Note that the gap value calculated within DFT is in general underestimating the true band 
gap value. 

In this paper we consider the sublattice symmetry breaking mechanism for a gap open- 
ing in a pristine doped graphene sheet and study the impact of gap upon some electronic 
properties of QPs. To investigate the influence of gap in the many-body properties of QP 
in graphene we use the random phase approximation (RPA) and the GqW approximation. 
It should be noted that a detailed analysis provided a framework for the microscopic eval- 
uation of the QP-QP interaction in the gapless graphene by means of the RPA was carried 
out by us in Ref. 49|] At the beginning, we review briefly the results of the ground state 



thermodynamic properties that we have already presented elsewhere. [50| Our new results 
are based on the QP self-energy properties in the presence of a gap opening in the electronic 
spectrum. From the self-energy we then obtain the QP energies, renormalized Fermi veloc- 
ity, spectral function which can be compared with ARPES spectra and finally the band gap 
renormalization of massive Dirac Fermions in doped graphene. We have shown that mass 
generating in graphene washes out a satellite band in the spectral function in agreement 



with recent experimental observations. 



43| 



This paper is organized as followed. In Section II we introduce our model Hamiltonian 
and then review some ground state properties of gapped graphene. In Section III we focus 
on the properties of imaginary and real parts of self- energy for gapped graphene and then 
calculate QP spectral function, renormalized Fermi velocity and band gap renormalization. 
Finally we conclude in Section IV. 



II. GROUND STATE THERMODYNAMIC PROPERTIES 

We consider the sublattice symmetry breaking mechanism in which the densities of par- 
ticles associated to on-site energy fia{b), for A(B) sublattice are different. The electronic 
structure of graphene can be reasonably good described using a rather simple tight-binding 
Hamiltonian, leading to analytical solutions for their energy dispersion and related eigen- 
states. The noninteracting tight binding Hamiltonian for vr band electrons is determined 
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by |38|, l39|, |40|, Q 



Ho = t ^(a|6i + c.c.) + /ia ^ aUi + f^b'^ h% 

i i i 

= t J2{alk + c.c.) + 5^(ala, - blk) + ^^(ala, + fetfe.) (1) 

i i i 

where the sums run over unit cells, t ~ 2.7 eV denotes the nearest neighbor hopping pa- 
rameter and ai{bi) is Fermi annihilation operator acts on A(B) sublattice. The second term 
in the noninteracting Hamiltonian breaks the inversion symmetry and causes to a band gap 
with value of 2A = \fia — fJ'b\ at the Dirac points. The last term is a constant and we 
left it out. The effective Hamiltonian at low excited energies lead to a 2D massive Dirac 
Hamiltonian, 7io = hvpa ■ k + Aas, where a are Pauli matrices and vp = 3ta/2h ~ 10^ m/s 
is the Fermi velocity where a ~ 1.42 A is the carbon-carbon distance in honeycomb lattice. 
The two eigenvalues of noninteracting Hamiltonian are given by = ±^ihvFky + for 
conduction band (+) and valance band (-) which is a fully occupied. In addition, the model 
Hamiltonian can be used as an approximated model for describing a graphene antidot lattice 



in the vicinity of a band gap with a small effective mass value 5^, or moreover used as an 



effective Hamiltonian for the intrinsic spin-orbit interaction in graphene where A = Aso 



is the strength of the spin-orbit interaction. [34, l35|, l36|, l37[ If fJ^a = l^b the Hamiltonian 



reduces to massless Dirac Hamiltonian with two chiral eigenstates having the conical band 
structures Ek = ±hvpk. 

We consider the long-range Coulomb electron-electron interaction. We left out the inter- 
valley scattering and use the two component Dirac Fermion model. Accordingly, the total 



interacting Hamiltonian in a continuum model at point is expressed as 57|, |58| 

^ = E ^L^O^k,. + ^.K^-q - (2) 

where \E'k^ = {ip'^ ^(k),ip^ ^(k)) is two component pseudospinors of the noninteracting 
Hamiltonian, S is the sample area, N is the total number operator and Vg = 2Tie^/eq is 
the bare Coulomb interaction where e is an average dielectric constant of the surrounding 
medium. The coupling constant in graphene is agr = gsOv^^ /^^vp where gs = = 2 being 
the spin and valley degeneracy, respectively. The coupling constant in graphene depends 
only on the substrate dielectric constant while in the conventional 2D electron systems 
is density dependent. The typical value of dimensionless coupling constant is 1 or 2 for 
graphene supported on a substrate such a SiC or Si02. 



A central quantity in the many-body techniques is the noninteracting dynamical polar- 
izability function ^(^^(q, zcj, /i) where /i is the chemical potential. The problem of linear 
density response is set up by considering a fluid described by the Hamiltonian, which 
is subject to an external potential. The external potential must be sufficiently weak for 
low-order perturbation theory to suffice. The induced density change has a linear relation 
to the external potential through the noninteracting dynamical polarizability function. This 
function is recently calculated along the imaginary frequency axis and it is given by jsol 



Ay 



Ay 



x„(sin sin ) 



(3) 

where a;± = y^l ± AA? /{e"^ + h?uj'^), y = ^/e^ + ti^ijj'^ and z{x) = (2x + ihuj)/eq. The Fermi 
energy of a 2D massive Dirac Fermion system is given by Ef = /i = a/ {hvpkpy + and the 



Fermi wavevector is = ^yAirn/ggg^ where n is the density of carriers. The noninteracting 
density of states (DOS) is determined by D{E) = gsgv\E\/27i}fvpQ{E'^ — A^) which is 
density dependent at the Fermi surface. It should be noticed that D{Ep) equals to m/27r/i^ 
in the conventional 2D electron gas system. Here, Q{x) is Heaviside step function. 

We now turn to present our first numerical results which are based on the noninteracting 
polarization function. The static polarization function as a function of wavevector for various 
gap values is shown in Fig. 1(a). The static polarization function in gapless case is a smooth 
function whereas a kink at g = 2kp occurs for gapped graphene and thus the derivatives of 
x'^'''*(q, 0, A 7^ 0) has a singular feature. The singular behavior is the source of several phe- 
nomena such as the Friedel oscillations and moreover the Ruderman-Kittel-Kasuya-Yoshida 
(RKKY) interaction which the later is absent in gapless graphene. In Figs. 1(b) and (c) we 
have plotted the dynamic polarization function as a function of frequency for wave vectors 
smaller and larger than q = 2kp, respectively. x*^'^^(q, zu;) tends to zero like uj~^ at large 
frequency region. 

The polarization function along the real u axis can be obtained by performing analytical 
continuation of Eq. ([3]). 50, l59|, l60|] In Fig. 2 we have presented the real and imaginary 
parts of the noninteracting polarization function as a function of frequency. Sharp cutoffs in 
the imaginary part of ^(^^(q, a;) are related to the rapid swing in the real part of x^^K^^'^)- 
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FIG. 1: (Color online) (a): The static noninteracting polarization function as a function of q for 
various A. The dimensionless noninteracting dynamic polarization at (b): q = O.bkp and (c): 
q = 2/cf as a function of uj for various A. 

These behaviors are in result of the fact that the real and imaginary parts of the polarization 
function are related through the Kramers- Kronig relations. Importantly, the sign change of 
the real part from negative to positive shows a sweep across the electron-hole continuum. At 
very large gap values, the polarization function of massive Dirac Fermions can be reduced to 
the polarization function (the Lindhard's function) of conventional two dimensional electron 
gas systems, as they are determined in Figs. 1 and 2. Consequently, we settle under situation 
that we can describe a range of band structures from the Dirac's cone (gapless graphene) 
to the parabolic (conventional semiconductors) band structure behavior by tuning the gap 
values from zero to a large value, respectively. We limited our calculations to the interme- 
diate values of of A and we thus expect wide range of the particular properties related to 
unique behavior of the polarization function. 



58 



We can calculate the total ground state energy of gapped graphene within RPA. [50|, 



The ground-state energies can be calculated using the coupling constant integration 
technique, which has the contributions Etot = £kin + ^xc- The kinetic energy per particle is 
given by €kin = 2(^| - A^)/3el. 



As discussed previously [50|, 158| we might subtract the vacuum energy contribution from 
the total energy, 6etot = ^totik-p) —Stotik-p = 0). Due to the number of states in the Brillouin 
zone must be conserved, we do need a ultraviolet cut-off kc, which is approximated by 
irk"^ ^ (27r)^/Ao, where Aq is the area of the unit cell. The dimensionless parameter A is 
defined as kjkp ~ {gsg^n'^VS/9mY/^ x lO^. 
In Fig. 3, we have shown the exchange- correlation energy in units of = hv-pkp, as a 
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function of n~^/^ in units of 10~^ cm for various A value. The exchange energy arises entirely 
from the antisymmetry of the many-body wave function under exchange of two electrons 
is positive while the correlation energy, the difference between the ground state energy and 
the sum of the kinetic energy and the exchange energy is negative. This has important 
implications on the thermodynamic properties can be calculated from the derivative of the 
ground state energy with respect to the density. The compressibility can be calculated from 
its definition, = n'^d'^ {n6etot) / dn? . Fig. 3(b) shows the ratio between the noninteractiong 
value, kq = 2/ne-p and the interaction value of compressibility as a function of The 
exchange tends to reduce the compressibility while correlations tends to enhance it. At 
large A, a minimum structure occurs at the inverse of compressibility behavior and we 
expect that at very large A, it starts at kq and reduces by increasing n"^/"^ behaves like the 
compressibility of the conversional 2D electron gas. 
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FIG. 2: (Color online) The gap dependence of the real and the imaginary part of the noninteracting 
polarization function as a function of u; for wavevectors (a), (c) q = 0.5/cf and (b), (d) q = 2.5kp 
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FIG. 3: (Color online) (a) Exchange-correlation energy and (b) compressibility as a function of 
n'^^l'^ ( in units of 10~® cm) for various A value. 

III. THE QP SELF-ENERGY AND THE SPECTRAL FUNCTION 

The generation of QPs in an electron liquid leads to two effects. First it induces a decay 
of a particle losing momentum via inelastic scattering which is determined by the imaginary 
part of self-energy and second is the renormalization of the dispersion relation of the carriers 
which is described by the real part of self-energy. 3fJeS'''^*(k, cj) is defined as the difference 
between the measured carrier energy hio^ and the energy of free particle, = S-Ek — E^. 
To satisfy causality, the real and imaginary parts of self-energy are related by a Hilbert 
transformation. In this section, we first derive the imaginary and the real part of QP 
self-energies and then calculate some important quantities such as a renormalized Fermi 
velocity, a spectral function and a band gap renormalization in the presence of a band gap 
value. These quantities are related to some important physical properties of both theoretical 
and practical applications like the band structure of ARPES, the energy dissipation rate of 



injected carriers and the width of the QP spectral function. 61 



62| 



In the GpW approximation, the self-energy of gapped graphene is given by (/? 
l/(A:ijT)) |63|: 

1 ^ /■ d^q 



S,,(k, 



s' 

+ 00 



(4) 



m=—oo 



where W{q,iQrn) = Vq/e{q,iflm) is the dynamical screened effective interaction and 
e(g, i^m) = 1 — ygX^^Hl^ i^m) is dynamical dielectric function in RPA. The overlap function 



for gapped graphene F**'(k, k + q) arises from the graphene band structure is given by [5 

^ -C/k-C/k+q 

It should be noted that F^"~^'(q = 0) = 0. However, in gapless graphene, intraband back- 
ward scattering should not be allowed, namely F^^^' {q = —2k, A = 0) = 0, as well as 
F'=-''{q = 0,A = 0) = 0. In Eq. (0]), Gf\k,iuj) = l/{iuj -^sk/h) is the noninteracting 
Green's function. Notice that in typical density of carriers in graphene namely n > lO^^cm"^, 
the Fermi temperature is about Tp = e-p/ks > lO'^K, and we therefore can eliminate temper- 
ature parameter in our calculations. To evaluate the zero-temperature retarded self-energy 
we perform the line-residue decompositions, E^'^*(k, c<j) = S^"^'^(k, cu) + S^'^'^(k, cu), where S'"^*^ 
is obtained by performing the analytic continuation before summing over the Matsubara 
frequencies, and T,'^'^^ is the correction which must be taken into account in the total self- 

n 

energy. [57| At zero temperature we have 

Si-(k, CO) = -J2 [ 7|^Ki^^^'(k, k + q) (6) 



X 



i2ny 

dn 1 1 

2ti e(q, i^l) u + - ^s' (k + q)/^ ' 



and 



X [e(c^ - ^,,(k + q)/h) - Qi-^Ak + q)/h)]. (7) 

The line contribution of the self-energy is purely real. The imaginary part of the self-energy 
has two contributions where 3mS'^*(k, cj) = Q'mS5^'^j^^j,jj(k, cj) + 53mZ15^''jj^^.pj.(k, u;), and real 
part of the self energy can be decomposed as 3?eS^]!:*(k, c^;) = S!j,'^'^(k, u) + 3?eS^^*\^^gj,(k, u) + 
5ReEWa(k,a;). 

For uj > and fixed q, the RPA decay process represents scattering of an electron from 
momentum k and energy uj to k -\- q and ^s'{k + q), with all energies in Eq. ([7]) measured 
from the Fermi energy of doped graphene. Since the Pauli exclusion principle requires that 
the final state is unoccupied, it must lie in the conduction band, i.e. s' = +1. Furthermore 
since the Fermi sea is initially in its ground state, the QP must lower its energy, i.e. ^s' < co, 
electrons decay by going down in energy. For u < 0, the self-energy expresses the decay 
of holes inside the Fermi sea, which scatter to a final state, by exciting the Fermi sea. 
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In this case the final state must be occupied so both band indices are allowed for s', and 
energy conservation requires that holes decay by moving up in energy. Since photoemission 
measures the properties of holes produced in the Fermi sea by photo ejection, only < is 
relevant for this experimental probe. 

In what follows, we calculate the intraband and interband contributions of self-energy. 
We have found the intraband term of residue part of self-energy as following for various 
values of the frequencies, 

SV%,,(k, u>0) = C dq dyUiy, q) 

^max(0, kp— k, k— /3) ^max(Ep, o_) 

pk+kp /■min(Ep, q+) 

SV^: i,,,,(k, A-Ej.<hLU<0) = -C dq dyUiy, q) 

^max(0, k— kp, /3— k) ^max(0, hw+Ej-p, o_) 
pk+kp /•min(Ep, 

EV%,(k, ^ < -(A + Ep)) = -C dq dyUiy, q) , (8) 

^max(0, k— kp) ^max(0, hw+'Ep, a.—) 

where 



e(q, u + EfTv) \/^k'^q^ - {y^ -El- q^f 
C = eV27reEfc, a± = ^/m4{k±qfT^ and p = ^h^u^ + h^v^4 + 2huEF. On the 
other hand, the interband contribution of residue part of the self energy is determined by 

r>fc+/3 ^min{o+,— fej— Ep) 

'max(0, k~/3) 

and eventually for the line contribution of self-energy we have 



rk+P /■min{o+,— fej— Ep) 

S-:,„,,,(k, nuj<-{A + Ef)) = -C dq dyf.{y, q) , (9) 

JmaxfO, k-fll >/o_ 



where 



S!r:Lra(k,^) = f^dq rrf0F++(q,q + k,A) P rffi^^^lM 

2 pkc r27T p+oo ((h Q n) 

= I I *F-(q.q + k.A) (10) 

and denotes an angle between k and q. Note that the real part of self-energy is kc 
dependent. 

Now we are in a situation that can calculate some important physical quantities. The QP 
lifetime or the single-particle relaxation time r, is obtained by setting the frequency to the on- 
shell energy in imaginary part of the self-energy, t^^ = Ts{k,C,sk/}T') = ||^5mE^*^*(k, ^^k/^) | 
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where rs(k, ^^k/^) is the quantum level broadening of the momentum eigenstate |sk >. This 
quantity is identical with the Fermi's golden rule expression for the sum of the scattering 



rate of a QP and quasihole at wavevector k. 57| From Eqs. 8 and 9, one can conclude that 



total contribution of the imaginary part of the retarded self-energy on the energy shell comes 
from the intraband term, 53mS5^*(k, ^k/^) = ^"^Sintra(^^> ^k/^)- j63|] In the case of gapless 



graphene, scattering rate is a smoot 



emission and interband processes. 



1 function because of the absence of both plasmon 
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65| However, with generating a gap and increasing 



the amount of it, plasmon emission cause discontinuities in the scattering time, similar to 



conventional 2D electron gas. 66|,l67| We have thus two mechanisms for scattering of the QPs. 



The excitation of electron-hole pairs which is dominant process at long wavelength regions 



63|, 



and the excitation of plasmon appears in a specific wave vector. As discussed previously 
in clean graphene sheets the inelastic mean free path reduces by increasing the gap whereas 
the mean free path is large enough in the range of the typical gap values 10-130 meV, and 
thus transport remains in the semi-ballistic regime. 

The many-body interactions in graphene as a function of doping can be observed by 
ARPES which plays as a central role to investigate QP properties such as group velocity 
and lifetime of carriers on the Fermi surface. ARPES is a useful complementary tool which 
capable of measuring the constant energy surfaces for all partially occupied states and the 
fully occupied band structure. The information of band dispersion and the Fermi surface can 
be elicited from those data measured in ARPES experiments. The relation of the Green's 
function to the single-particle excitation spectrum in the interacting fluid is expressed by 



its spectral function. 



following expression 57 1 



'he spectral function is related to the retarded self-energy by the 



where 5Sg'^*(k, cj) = S^°*(k, cj) — T,l'^^{kp,0), and then ARPES intensity can be described by 
/(k, cij) = A{k, u!)n{u!) , where n{u) is the Fermi-Dirac distribution. The spectral function is 
the Lorentzian function where 3?eS specifying the location of the peak of the distribution, 
and |53mS| is the linewidth. The amplitude of the the Lorentzian function is proportional 
to l/|$5mS|. This quantity is the distribution of energies hw, in the system when a QP 
with momentum k, is added or removed from that. For the noninteracting system we get 
A^'^'' (k, u) = 6{uj — C,(k)/h). The Fermi liquid theory applies only when the spectral function 
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at the Fermi momentum A^^\k = kp,Lj), behaves as a delta function, and has a broadened 
peak indicating damped QPs at k ^ kp. 

To progress of the interband single particle excitation and plasmon effects on the Q'mE^, 



we must study the retarded self-energy on the off-shell frequency which is u; 7^ ■Csk/^- 491.168 1 
This quantity gives the scattering rate of a QP with momentum k and kinetic energy hw+Ep. 
The scattering rate or the linewidth raising from electron-electron interactions is anisotropic 
and varies significantly via wavevector at a constant energy. The imaginary part of self 
energy shows the width of the QP spectral function. 

In Fig. 4 we have shown the absolute value of the imaginary part of the self energy in 
unit of Ep for various gap values. It would be noticed that there is an area of frequency 
is associated to the gap value, 2A in which no QP could enter in. In this case, there is 
a gap in the 53mS between ^-,fc=o and ^+^k=o- We see that 53mS+ vanishes as u"^ for lj 
tends to zero, a universal properties of normal Fermi liquid. Moreover, at large frequency, 
Q'mS_|_ tend towards to u linearly. Except from the Dirac point, the conduction band 
'ismTj^ peaks broaden because of the dependence on scattering angle of ^(k + q). For low 
energy, only intraband single particle excitation contributes to QmT, up to Ep and then the 
interband single particle excitation contribution increases sharply about Ep. The interband 
contribution increases with increasing the gap values. 

To evaluate the scattering rate in interband channel, we have shown QmHinterimtra) as 
a function of frequency in Fig. 5. The intraband contribution of the imaginary part of 
self energy associated to scattering rate of QP in the intraband contribution increases with 
increasing the gap values while the interband contribution reduces, as we physically expected. 
Moreover, by increasing of the electrons in the conduction band the interband scattering rate 
reduces whereas the intraband scattering contribution increases. The gap value suppresses 
the scattering rate at = —ep. 

In Fig. 6 we have plotted the real part of self-energy in unit of ep as a function of the 
energy for various gap values. Notice again that the real part of residue self-energy has a 
gap which is associated the feature calculated in the imaginary part of self-energy. The line 
part of self energy is a continues curve and then we have a jump near to the boundary of gap 
values in the 3?eS for gapped graphene. A kink around Ep is associated to the interband 
plasmon contribution and it is broaden due to the gap value. This feature affects noticeably 
in the interacting electron density of states. 
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FIG. 4: (Color online) The absolute value of the imaginary part of retarded self-energy ( + channel) 
as a function of (a) uj and (b) k, for the various energy gaps. 
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FIG. 5: (Color online) (a) The intraband and (b) the interband contributions of the imaginary 
part of self-energy ( + channel) as a function of uj for the various energy gaps at k = 0.25kF. 



As discussed before 



49 



68l | in a zero temperature and disorder free gapless graphene, the 
peaks of the spectral function correspond to the nearly solutions of Dyson's equation in which 
the quasiparticle excitation energies are obtained by = ^+ + 3?e5S^^*. The intersection of 
3fJeE and the lines E — C,+ indicates a satellite long wavelength plasmaron peak related to the 
electron-plasmon excitation due to the long-range electron-electron Coulomb interaction and 
the Dyson equation with 53mS = corresponds to a QP peak related to the single particle 
excitation. Importantly, in the presence of gap values, the plasmaron peak suppressed. 
In Figs. 7(a) and (b) we have shown the energy distribution curves (EDC) and momentum 
distribution curves (MDC), respectively. In the presence of gap values, as shown in Fig. 7(a) 
there is only the single QP peak. 

The valance band self-energy contributions are shown in Fig. 8. There is an area of 
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FIG. 6: (Color online) The real part of retarded self-energy ( + channel) as a function of (a) lo 
and (b) k, for the various energy gaps at A = 100. The 3^eS_|_ are measured from the interaction 
contribution of the chemical potential 3fteS+(/cF,u; = 0). 
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FIG. 7: (Color online) The QP spectral function for + channel as a function of u; (a) and k (b), 
for the various energy gaps at A = 100. 

frequencies is associated to the gap value in which no QP could exist in 53mS_ exactly 
the same as the conduction band. The s = +1 and —1 peaks in SJmSs in Figs. 4 and 8 
separate at finite k because of chirality factors which emphasize k and q in nearly parallel 
directions for conduction band and k and q in nearly opposite directions for valance band 
states. Consequently, at finite A, the QP peak of A_{k^iS) which is broaden shifts toward 
the left in the opposite behavior of Aj^^k^iS). These feature have significant effects in the 
interacting electron density of states. 



It is essential to note that the satellite band which is theoretically predicted [49| for gapless 
garphene has not been seen in experiments. There are several reasons that could wash out 
this feature. For example, the plasmon damping, disorder effects, electron interactions with 
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function as a function of to (c), for the various energy gaps at A = 100. 



> 




FIG. 9: (Color online) The renormalized velocity (+ channel) as a function of the energy gap for 



various densities at agr = 1. 



the buffer layer and importantly the effect of gap at Dirac point. 

One of the important information which can be extracted from ARPES spectra is the 
renormalized Fermi velocity v* . A consequence of the interaction is a Fermi velocity renor- 
malization from the backflow of the fluid around a moving particle. The density of states at 
the Fermi energy is also changed. The QP energy measured from the chemical potential of 
interacting system 5e^^ , can be calculated by solving self consistently the Dyson equation 



571 ] In the isotropic systems the QP energy, depends on 



the magnitude of k. Expanding 6e^^ to first order mk — k-p we can write 5e^^ ^ hv*{k — kp) 
which effectively defines the renormalized velocity as hv* = d5e^^ / dk\k=kY- From the Dyson 
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equation we can calculated the renormalized Fermi velocity as [63|, |69|, |7 

^ BfE^' + {hvs)~'dk^e[6^T{k,uj)]\^=o,k=k^ 



v.. 



where f , 



1 - n-i(9^3?e[(5S^/t(k,w 
svp. It is found before 



\uj=0,k=kp 



(13) 



631, 



69 



7G 



71 



72] that electron-electron interaction 



increases the renormalized Fermi velocity in gapless graphene sheets which this behavior is 



in contrast to conventional 2DES. 



66 



73| 



Fig. 9 shows the renormalized Fermi velocity in unit of the bare Fermi velocity as a 
function of band gap for various carrier densities. The renormalized Fermi velocity decreasing 
with increasing the gap value, v* is density independent after A = O.Sep which is in good 



24 



25| 



76 



77| The BGR 



agreement with recent experiment observation. ^ ^ 

Finally we calculated a band gap renormalization (BGR). 
for conductance band is given by the QP self-energy at the band edge, namely BGR = 
3fJeS!^'(k = 0,uj = {A — EF)/h). Fig. 10 shows the BGR for the various gap values as a 
function of the electron density. The BGR decreases by increasing of the electron density 
and in the small energy gap values, it is less density dependent respect to large energy 
gap values. In gapless case, we have obtained a induced band gap or kink due to many- 



body electron-e 



density. [49 



68 



ectron interactions and it tends to a constant with increasing the electron 



71 



DFT calculation. [Til 



73] This feature is in agreement with the results obtained within ab intio 
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FIG. 10: (Color online) The band gap renormalization as a function of electron density for various 
energy gap at Ugr = 1. 
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IV. SUMMERY AND CONCLUSION 



We have revisited the problem of the microscopic calculation of the QP self-energy and 
many-body effective velocity suppression in a gapped graphene when the conduction band 
is partially occupied. We have performed a systematic study is based on the many-body 
GqW approach that is established upon the random-phase-approximation and on graphene's 
massive Dirac equation continuum model. We have carried out extensive calculations of both 
the real and the imaginary part of the QP self-energy and discussed about the interband 
and intraband contributions in the scattering process in the presence of gap value. We have 
also presented results for the effective velocity and for the band gap renormalization over a 
wide range of coupling strength. Accordingly, we have critically examined the merits of the 
gap values in dynamical QP properties. 

Most feature of mass generating in graphene is the washing out of the plasmaron peak in 
the spectral weight. Increasing of the gap value makes density independent behavior of the 
renormalized Fermi velocity. We have shown that the band gap renormalization in gapped 
graphene decreases by increasing the carrier density at large A. This is in contrast with 
the gapless case in which many body electron-plasmon interactions induce a very small gap 
in band structure. These distinct features of the massive Dirac's Fermions are related to 
mixing of the chiralities and reduce of the interband transitions in graphene sheets. 
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